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Abstract. While ion heating by clastic electron- ion collisions may be neglected 
for a description of the evolution of freely expanding ultracold neutral plasmas, the 
situation is different in scenarios where the ions are laser-cooled during the system 
evolution. We show that electron-ion collisions in laser-cooled plasmas influence the 
ionic temperature, decreasing the degree of correlation obtainable in such systems. 
However, taking into account the collisions increases the ion temperature much less 
than what would be estimated based on static plasma clouds neglecting the plasma 
expansion. The latter leads to both adiabatic cooling of the ions as well as, more 
importantly, a rapid decrease of the collisional heating rate. 
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Recently, the field of ultracold (T <C IK) neutral plasmas and Rydberg gases has 
attracted attention both experimentally [UI21E1 and theoretically [U EJUS 0011]]. One of 
the motivations of the experiments [TJ Ell ^ s the creation of a so-called strongly coupled 
plasma, where the Coulomb coupling parameter r ; = e 2 /(ak^Tj) is much larger than 
unity (where a = (47rpi/3) -1 / 3 is the Wigner-Seitz radius and Tj is the temperature of 
the ions). In such a case, interesting ordering effects such as Coulomb crystallization 
into cubic or shell-structure lattices can be observed. In the experiments strongly 
coupled ions can not be observed since the temperature of the initially cold ions rises 
quickly on the timescale of the inverse ionic plasma frequency u~j = \fm\j [Ame 2 p\) 
due to disorder- induced heating jH El Ej . Different ways to overcome this heating effect 
have been proposed jH HTJ O EH EH IH3- In [T21IT3] . e.g., it has been suggested that 
continuous laser cooling during the plasma expansion may considerably increase the 
achievable Coulomb coupling parameter of the ions. 

Using a hybrid molecular dynamics approach (EH , we have predicted that Coulomb 
crystallization in an unconfined ultracold neutral plasma can indeed be obtained if the 
plasma ions are laser-cooled during the evolution of the system |14j . In this theoretical 
approach, the electronic component of the plasma is treated as a fluid while ions and 
recombined atoms are described on a full molecular dynamics level. More precisely, an 
adiabatic approximation is made for the electrons, which are assumed to be distributed 
according to the mean-field potential generated by ions and electrons, calculated in a 
self-consistent way. Their velocity distribution is assumed to be of Michie-King-type 
[T7] which takes into account deviations of the quasi-equilibrium state from a Maxwell- 
Boltzmann distribution due to the finite depth of the potential well generated by the 
ions. The latter leads to an evaporation of a fraction of the electrons in the initial stage 
immediately after the plasma formation, which is taken into account using the results 
of pQ . The ions, on the other hand, are propagated under the influence of all other ions 
and the electronic mean-field. Finally, inelastic collisions, i.e. three-body recombination, 
electron-impact ionization and electron-impact (de)excitation, are included on the basis 
of a Monte Carlo treatment. Elastic electron-ion collisions, however, are neglected in 
this approach, as in other approaches for the description of ultracold plasmas [TH] . 

This neglect of elastic electron- ion collisions in the dynamics is well justified for 
the description of freely expanding ultracold plasmas as created in [I] because of a clear 
separation of timescales. Due to the large mass ratio between ions and electrons, the time 
necessary for equilibration of the electron and ion temperature is typically of the order of 
a few milliseconds, while the plasma expansion takes place on a microsecond timescale. 
Hence, elastic electron-ion collisions only increase the ion temperature by some milli- 
Kelvin during the experimental observation time £Q . This amount of heating is negligible 
compared to the initial temperature increase due to disorder-induced heating, which 
rises the ion temperature to about one Kelvin. Hence, the time evolution of the ion 
temperature which enters T; and determines the degree of coupling of the ions is mainly 
determined by this disorder-induced heating. 

However, this situation changes in the scenario of [T3], where additional laser cooling 
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of the ions compensates the disorder- induced heating and keeps the ionic temperature on 
a milli-Kelvin level. In this case, the ion temperature is not only increased by electron- 
ion collisions, but is also driven to the Doppler temperature, i.e. the limiting temperature 
T c for laser cooling, due to the coupling to the radiation field. The final ion temperature 
is now determined by the balance between the collisional heating rate and the laser- 
cooling rate ^21 EI]- Since T c is typically of the order of one milli-Kelvin, i.e. of the 
same order of magnitude as the amount of collisional heating, the ion temperature may 
considerably increase through electron-ion collisions. Thus, electron-ion collisions could 
significantly decrease the Coulomb coupling parameter achievable by laser-cooling of 
expanding ultracold neutral plasmas. In the following, the influence of such collisions on 
the onset of Coulomb crystallization during the plasma expansion is critically reassessed. 

The system under study is the same as that of ^1] , namely an unconfined, ultracold 
neutral plasma under the additional influence of a cooling laser. The plasma is treated 
on the basis of the hybrid-MD method outlined above and described in detail in [15] . 
Laser cooling is modelled by adding a Langevin force, F coo i = — m^Vi + y/2f3k-BT c m£, 
to the ion equation of motion, where Vi is the ion velocity, £ is a stochastic variable with 
(£) = 0, (£(£)£(£ + t)) = 3<5(t) and the cooling rate (3 and the corresponding Doppler 
temperature T c are determined by the properties of the cooling laser |19j . The weak 
electron coupling, which justifies our hybrid treatment of the plasma dynamics, also 
allows us to describe electron-ion collisions with a Boltzmann-type collision operator, as 
we shall discuss below. The implementation of electron-ion collisions is very similiar to 
the treatment of electron-electron collisions used in ^H] and has been introduced before 
in [201 • For an i° n a t position rj with velocity v ; , the rate of electron-ion collisions is 

K eic = [ dv e dB sin6 ( |vj - v e | /e(r i( p e ) . (1) 



Hence, P e j C = At K eic is the probability that a collision will occur during a timestep 
At. Numerically, the integration over the electronic velocity distribution in equation 
is done via a Monte Carlo procedure. At each timestep, an electronic velocity v e is 
chosen randomly according to a Maxwell-Boltzmann distribution with temperature T e . 
The probability for an electron-ion collision is then given by 

r 1 ( da\ ne 4 
P cic (v c ) = At p c ( ri K / — rf(cos6) = At Pc (r i )v e ^ 1 (A - 1) , (2) 

J -i yafcJ j eff m e V e 

with In A = In (^/3/T^ 2 ^j the so-called Coulomb logarithm and 

—) 9 > 6 ■ 

deJconi u - Umin (3) 

o e < e r 




'mm 



where (^§) c } is the Rutherford cross section for Coulomb scattering [21] and 6 m i, 
2arcsinA -1 / 2 [22]- In the derivation of equation ([2]), the approximation m e /m; 
0, vi — > has been made (see below). 
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Figure 1. Time evolution of the ionic Coulomb coupling parameter with (solid) and 
without (dashed) the inclusion of ion heating by electron-ion collisions for a plasma 
of 20000 Be ions with p ia = 2.3 x 10 8 cm" 3 , r c0 = 0.05, f3 = 0.10 u p ,i(t = 0) and 
T c = lmK. 

In the Monte Carlo treatment outlined above, an electron-ion collision takes place 
with a probability P C i C (v ). In this scattering; angle G is determined according; 

to the probability measure P(G) oc (4^1 , leading to the prescription 



where ( is a random number distributed uniformly in [0, 1]. With this, the change of 
ionic momentum is given by Ap; = m e v c (1 — cos Q) j, where j is a vector denoting 
the random direction of the momentum transfer chosen uniformly on the unit sphere. 
Finally, energy conservation is restored by adjusting the electronic temperature T e . 

As mentioned above, the derivation of equation (0) implies the limit m e /mi — > 0. 
In this case, individual collision events satisfy momentum conservation, but energy 
conservation is violated and has to be corrected manually by adjusting the electron 
temperature. The error introduced by this approximation can be shown to be of the 
order of the temperature ratio Ti/T e [23! • Hence, as long as the ionic temperature is small 
compared to the electronic temperature, which is the case over the whole time of the 
experiments under consideration, the corresponding corrections can be neglected. The 
method described above thus allows for an efficient simulation of electron-ion collisions 
without the need for time-consuming transformations between laboratory and center- 
of-mass coordinate systems. 

In order to check the influence of elastic electron-ion collisions on the plasma 
dynamics, we have simulated the expansion of a plasma with an initial electronic 
Coulomb coupling parameter of T e {t = 0) = r e o = 0.05, consisting of 20000 Be 
ions at a density of = 0) = p i0 = 2.3 x 10 8 cm" 3 , cooled with a damping 
rate of (3 = 0.10 uj p ^(t = 0) and a Doppler temperature of 1 mK. As can be seen 
in figure ^ the inclusion of elastic electron-ion collisions significantly reduces the 
ionic coupling parameter which can be achieved in this case. However, the value Tj 
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Figure 2. Distribution of scaled interionic distances after a time of t — 40 /is, 
compared to the calculated pair-correlation function of an OCP at Ti = 400. The 
initial-state parameters are the same as in figure ^ 



obtained after several lis is still considerably larger than the critical value r\ c pa 174 
[21] for crystallization. The good agreement found between the distribution of interionic 
distances, calculated as described in [Hj, and the pair-correlation function of a one- 
component plasma at r ; = 400 [2H| also confirms that the system evolves well into the 
strongly coupled regime (figure GJ). We note here that Debye shielding of the ion- ion 
interaction, which is not included in the present considerations, tends to increase the 
crystallization limit. However, as shown in 26J, even for a comparably large electron 
coupling of T e = 0.2 the crystallization limit is increased to pa 198 only, hence the 
inclusion of Debye screening does not significantly alter the results obtained with the 
present model. 

In fact, the simulation results show that the final I 1 ! is much larger than what 
one would expect from a static estimate by equating the heating 7 e i C T'e resulting from 
electron- ion collisions to the cooling 2/3 (7] — T c ) of the laser cooling as done in [T^ fTT?j . 
From the procedure described above, 7 e i C is obtained as$ 

7eicT e = 7J- J dv e dO sinG^y^ (1 - cos6) 2 (^^j ^e/e(i"i, Pe) , (5) 

leading to the Landau-Spitzer expression for the average heating rate |2Zj 

where uJ Vfi is the electronic plasma frequency. The validity and extensions of the Landau- 
Spitzer formula, which was originally derived for weakly coupled plasmas, have been 
discussed in several publications [23 EE EE EDI • In |2Z1 , ion heating by binary collisions 
has been studied without employing the cutoff-procedure equation (JSJ). For T e < 0.25, 

| Again, this implies the limit Ti/T c — > which is well fulfilled over the whole timescale of the 
experiments as described above. 
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Figure 3. Time evolution of the ionic temperature for a plasma of 100000 ions with 
p i0 = 2 x 10 7 cm~ 3 , r o0 = 0.08 and (3 = 0.15 w p ,i(f = 0), using different levels of 
approximation (see text). The arrow shows the static estimate equation Q for the 
final ion temperature, the dotted line in the inset equation 10. 



their numerical results are reproduced by equation © to within 10%. Moreover, 
collective effects, which largely decrease the relaxation rate at strong coupling [2E], were 
found to be negligible even for arbitrarily strong ion-ion coupling, as long as the electron- 
ion coupling is weak and the resonance of the ion excitation spectrum lies far below that 
of the electrons [30! . Therefore, under the present conditions, which correspond to the 
parameter regime studied in [3U|, the Monte Carlo treatment described above yields an 
adequate description of the electron-ion temperature relaxation process. 

Balancing collisional heating with the laser cooling one might expect a quasi- 
equilibrium state with a final ion temperature of 



in the optimal limit that T c — > 0. For the parameters used in figure ^ equation (|5|) 
predicts a final T; « 50, which is much smaller than what is observed in the simulation. 

In order to trace the origin of this discrepancy, we have performed a second set of 
simulations. In one simulation, we have set the Doppler temperature T c equal to zero 
and have turned off the initial correlation-induced heating of the ions by propagating 
them in the framework of a particle- in-cell method JB], i>e. on a mean-field level, rather 
than by a full MD simulation. By doing this, only the competition between the laser 
cooling, the heating by electron-ion collisions, and the adiabatic expansion determines 
the evolution of the ion temperature. In the second type of simulation, T c is also set 
to zero and all particle interactions except the binary electron-ion collisions causing the 
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leading to a Coulomb coupling parameter 
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ion heating are neglected, which results in a basically static ionic density§. 

The resulting time evolution of the ionic temperature for a plasma of 100000 ions 
with p i0 = 2 x 10 7 cm~ 3 , r e0 = 0.08 and (3 = 0.15 u P) i(t = 0) is shown in figure El In the 
static case (dashed line in figure EJ), the temperature shows an initial linear rise with 
a slope corresponding to the collisional heating rate 7 e i C T'e and then quickly saturates 
to a value which is well described by equation (JJJ) (marked by the arrow in the figure). 
However, in the case of an expanding plasma (solid line in figure El) the temperature 
is drastically reduced. While the initial rise of the ion temperature stays the same, 
the temperature starts to drop down already at a relatively early time, which at this 
stage is mainly caused by adiabatic cooling of the ions due to the plasma expansion. At 
later times, the ion temperature is driven to its steady-state value given by the balance 
between the collisional heating and the laser cooling (equation (JZJ)). The long-time 
behaviour of the ionic temperature can be estimated from the plasma dynamics discussed 
in where it was found that laser cooling strongly alters the expansion behaviour, 

leading to an increase of the plasma width a according to a oc tV*. It was also shown 
there that the adiabatic law for the selfsimilar plasma expansion, a 2 T e = const., still 
holds, which gives together with equation for T c — > 



while the ionic Coulomb coupling parameter increases according to T; oc l/(o~T e ) oc t 1//4 . 
As can be seen in the inset of figure El equation (jOJ is well reproduced by the numerical 
simulation. Hence, it is indeed this adiabatic cooling of the expanding plasma together 
with a decreasing collisional heating rate 7 e icT e which leads to much lower temperatures 
and consequently to much larger Coulomb coupling parameters than those predicted by 
a static estimate. 

The simulations described above demonstrate that electron- ion collisions do not 
prohibit achieving ionic Coulomb coupling parameters above the crystallization limit. 
However, for a given set of initial-state parameters, Fj can be significantly reduced by 
collisional heating during the plasma evolution, leading to an additional restriction of 
the initial-state parameters for which strongly coupled states can be achieved. Since 
it was found that the degree of spatial order in the system sensitively depends on 
initial conditions such as number of ions, electron temperature etc., one may wonder 
whether the long-range ordering into shell structures described in ^3] is still observable 
in simulations taking into account the electron-ion collisions. This is demonstrated in 
figure Elfor a plasma of 50000 ions with a density of p;o = 2 x 10 7 cm -3 and T e o = 0.08, 
cooled with a damping rate of 0.17ui Pt i(t = 0) and a limiting temperature of T c = 0.8mK. 
The ordering into concentric shells is clearly visible in the radial density after a time of 
t = 175 fis. 

In summary, we have presented simulations for laser-cooled ultracold neutral 

§ The slow diffusion of the ion distribution due to the finite ion temperature is negligible during the 
time of the simulation. 
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Figure 4. Radial density of a plasma of 50000 ions with initial-state parameters 
p i0 = 2 x 10 7 cm- 3 , r c0 = 0.08, (3 = 0.17w P)i (i = 0) and T c = 0.8 mK after t = 175ms, 
demonstrating the radial ordering in the inner plasma region. 



plasmas including elastic electron-ion collisions. It is found that these collisions can 
have a significant influence on the final ionic temperature, and hence the degree of 
correlation, achievable under given experimental conditions. Yet, with a proper choice 
of initial conditions rather large Coulomb coupling parameters can still be achieved. 
As we have shown, this is due to the fact that the expansion of the plasma drastically 
reduces the ion temperature, due to adiabatic cooling of the ions as well as a decreasing 
heating rate during the course of the plasma evolution. Therefore, the resulting Coulomb 
coupling parameters are much larger than one would expect from static estimates as in 
[121 HBj- We thus conclude that our predictions made in [T3] about the possibility to 
observe Coulomb crystallization in such a system remain valid. In particular, we have 
shown here that the electron-ion collisions do not preclude the build-up of long-range 
order and the formation of the shell structures predicted in The experimental 

realizability of the discussed scheme has been demonstrated recently in [TO]- There, the 
ion temperature of an ultracold Strontium plasma has been measured by driving the core 
transition of the plasma ions during the gas expansion. Since the same transition which 
allowed for Doppler imaging can be used to cool the ions, these experiments provide the 
first step towards a realization of strongly coupled ultracold plasmas in the laboratory. 

We would like to thank T.C. Killian and F. Robicheaux for valuable discussions 
and for drawing our attention to the importance of electron-ion collisions. Financial 
support from the DFG through grant R01157/4 is gratefully acknowledged. 
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